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We construct an exactly solvable spin-orbital model on a decorated square lattice that realizes an SU(2)- 
invariant Majorana spin liquid with parton Fermi surfaces, of the kind discussed recently by Biswas et al.[ Phys. 
Rev. B. 83, 245131(201 1)]. We find power-law spin correlations as well as power-law spin-nematic correlations 
with the same dominant l/|r| 3 envelope. The model is solvable also in the presence of Zeeman magnetic field. 
One fermion species carries S z = quantum number and its Fermi surface is not altered in the field, while 
the Fermi surfaces of the other species evolve and can disappear. In particular, we find an interesting half- 
magnetization plateau phase in which spin excitations are gapful while there remain spinless gapless excitations 
that still produce metal-like thermal properties. In the fully-magnetized phase, the model reduces to the one 
proposed by Baskaran et al, | arXiv:0908. 1614 ^3] in terms of the orbital degrees of freedom. 



I. INTRODUCTION 

Gapless quantum spin liquids (QSL)i^& are perhaps some 
of the most intriguing fractionalized phases. Much in- 
terest in these is motivated by recent experimental re- 
alizations in two-dimensional (2D) organic compounds 
EtMe 3 Sb[Pd(dmit) 2 ]2 and K-(ET) 2 Cu 2 (CN) 3 .i^ One pro- 
posal with Gutzwiller-projected Fermi sea wave functio n 29 ' 30 
is an appealing candidate but does not appear to be able to cap- 
ture all experimental phenomenology. Searching for alterna- 
tives, many possible proposals have been presented^ 3 - 4 . Very 
recently, Biswas et ali^i proposed an SU(2)-invariant Majo- 
rana QSL, which we find fascinating and in need of more at- 
tention. Motivated by this proposal, here we want to realize 
such long-wavelength QSL in an exactly solvable microscopic 
model. Following the route discovered by Kitaev^ and gener- 
alized to produce many other exactly solvable models.^— in 
particular with SU(2) spin invarianc o 52 ' 54 or with parton Fermi 
surfaces ^ 47 ' 49,53 we find a Kitaev-type model with both SU(2)- 
invariance and parton Fermi surfaces. 

Our model is realized using both spin- 1/2 and orbital de- 
grees of freedom^iSi at each site of a decorated square 
lattice. 47 . The system can be reduced to three species of free 
Majorana fermions coupled to background Z 2 gauge fields 
such that it is exactly solvable and parton Fermi surfaces 
are realized. We formulate long wavelength description in 
terms of an occupied Fermi pocket of three complex fermions 
(f x > f v > f z ) that transform as a vector under spin rotation. 
For general illustration (and also for preventing possible pair- 
ing instabilities away from the exactly solvable limit), we con- 
sider a model that lacks time-reversal and lattice inversion 
symmetries. Because of the exact solvability, we can learn 
much reliable physics information about such Majorana QSL. 

Specifically, we study spin correlations and spin-nematic 
correlations in our model. The main result is that these 
correlations have the same dominant power-law behaviors 
with 1 / 1 r | 3 envelope in real space and oscillations at incom- 
mensurate wavevectors which form what we call singular 
surfaces 1 ' 3,8 in the momentum space. Because of the Z 2 nature 
of the QSL and the absence of the time-reversal and inversion 
symmetries, there are additional non-trivial ±(kFR. + kpz,) 
and ±2k^ critical surfaces besides the more familiar km — 
\s.fl surface in the correlations. 



The model is still exactly solvable in the presence of Zee- 
man magnetic field. An interesting property is that the Zee- 
man field only couples to the f x and f y fermions while the 
f z fermion remains unaltered and therefore the f z Fermi sur- 
face remains and always gives gapless excitations. We cal- 
culate the magnetization as a function of magnetic field. In- 
terestingly, there is a plateau phase in which the spins are 
half-polarized with short-ranged spin correlations while the 
Fermi surface of f z still exists and gives gapless excitations, 
which can be detected using local energy operator like bond 
energy^ 

The paper is organized as follows. In Sec.[[l] we define the 
model on the decorated square lattice and solve it and discuss 
qualitative properties of the spin liquid phase. In Sec. [HI] we 
define the spin correlation functions and spin-nematic correla- 
tion functions. In Sec. lIII Al we provide a theoretical approach 
to describe the long-distance behavior of the correlations. In 
Sec. lIII B[ we present exact numerical calculations of the spin 
correlations and spin-nematic correlations. In Sec. |IV] we 
consider our model in the presence of the Zeeman magnetic 
field and specifically calculate the magnetization curve as a 
function of the field. We conclude with some discussion. 



II. SU(2)-INVARIANT MAJORANA SPIN LIQUID WITH 
STABLE FERMI SURFACES 



Motivated by the ideas from Baskaran et al£L, Yao et al^, 
and WangS, we construct an exactly solvable Kitaev-type 
model including both orbital and spin degrees of freedom with 
spin-rotation invariance. The Hamiltonian is 



H = H + H TRB + Ko Wo + Kq Wq , (1) 

o 



where, 
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The graphical representation of the model is shown in Fig. Q] 
At each site of the decorated square lattice, there are spin 
and orbital degrees of freedom. Ho is a Kugel-Khomskii-like 
Hamiltonian with a being the spin- 1/2 Pauli matrices and f 
being the Pauli matrices acting on the orbital states^Si The 
site labels in Eqs. (O-© are shown in Fig. Q] Htrb rep- 
resents an additional Time-Reversal-Breaking (TRB) interac- 
tion in the small diamonds 5 - (in principle, all four terms in the 
square brackets can have independent couplings). The reason 
for introducing the TRB and allowing different J x couplings 
in Ho that break the lattice point group symmetries is to avoid 
worrying about Cooper pair instabilities of the parton Fermi 
surface away from the exactly solvable limit. 

In addition, there are two types of elementary plaquettes 
(square and octagon) in the decorated square lattice (Fig. [T}, 
and two types of local conserved operators, Wo for the 
squares and Wq for the octagons in Eq. (0]i and Eq. (|5). The 
plaquette operators W p commute among themselves and with 
all other terms in the Hamiltonian and the K p terms are added 
to stabilize particular flux sector (see Fig.[T}. 

Introducing Majorana representation of spin-l/2^r— we 
write the spin and orbital operators as 



(8,7 



0,7 



(6) 
(7) 



On each site j of the decorated square lattice, we realize 
the physical four-dimensional Hilbert space using six Ma- 



jorana fermions c x ,. 



straint _D 7 



y 



-ic 



<l ' r <l ! '. and (/ 



with the con- 
1 (namely, for any physical 



state |$) p hys, we require Dj\<f>} pbys = |$) p h y s)- Therefore, 



phys = icfd?\$) phya 
In terms of the Majoranas, the Hamiltonian can be 




FIG. 1 . Graphical representation of the exactly solvable Kitaev-type 
model and its solution in the zero flux sector. The c x ' v,z Majoranas 
propagate with pure imaginary hopping amplitudes specified by the 
couplings J x , J y , J z , J x , J v , and J z ; the signs in our chosen 
gauge are indicated by the arrows. 



rephrased as 
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where Ujk = —idjcfy for A-link (jk). Following familiar 
analysis in Kitaev-type models, we observe that in the en- 
larged Hilbert space, iijk commute among themselves and 
with the Hamiltonian, and we can proceed by replacing them 
by their eigenvalues ±1 and interpreting as static Z 2 gauge 
fields. The W p terms, with K p > and assumed to be suf- 
ficiently large, can be used to stabilize the sector with zero 
fluxes through all elementary plackets, and this can produce 
parton Fermi surfaces.-^ In our work, we fix the gauge by tak- 
ing Ujk = 1 for bonds j — > k as shown by the arrows in Fig. Q] 
There are four physical sites per unit cell, so for each species 
c a , a = x,y, z, there are four Majoranas per unit cell. From 
now on, we replace the site labeling j with j = {r, a}, where 
r runs over the Bravais lattice of unit cells of the decorated 
square network and a runs over the four sites in the unit cell. 
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The Hamiltonian can be written in a concise form, 

= C r,a-Ar,a;r' ,a'C r i a , . (12) 

a ((r,o),(r',o')> 

There is translational symmetry between different unit cells, 
and Ar,a;r',a' = Aa' (r - r'). 

In order to give a concise long-wavelength description, it 
will be convenient to use familiar complex fermion fields. To 
this end, we can proceed as follows. For a general Majorana 
problem specified by an antisymmetric pure imaginary ma- 
trix Ajk, we diagonalize Ajk for spectra, but only half of 
the bands are needed while the rest of the bands can be ob- 
tained by a specific relation and are redundant. Explicitly, for 
a system with 2m bands, we can divide them into two groups. 
The first group contains bands from 1 to m with eigenvector- 
eigenenergy pairs {vb,k, e&,fe}> where b = 1, 2, . . . , m are 
band indices, and the second group contains bands from m+1 
to 2m related to the first group, {vb>= m +b,k, eh'=m+fe,fe} = 
{vt _fc 7 — e&.-fc}- Using only the bands with b = 1 to m, we 
can write the original Majoranas in terms of usual complex 
fermions as 



c a (r,a) = \/— J2 E [e Jkr ^,k(a)/ b Q (k) +H.c] ,(13) 



6=1 keB.Z. 



where N uc is the number of unit cells, and the complex 
fermion field / satisfies the usual anticommutation relation, 
{/ b at (k),/«'(k')} = <W<W<W. Note that in this SU(2)- 
invariant model, the eigenvectors for each spin species are the 
same, k = Ub.k- In terms of the complex fermion fields, the 
Hamiltonian becomes 



n = J2 E 2e ^ k ) 



6=1 keB.Z. 



/ b at (k)/6 Q (k) 



(14) 



In the present case, 2m = 4 and therefore two bands are 
sufficient to give us the full solution of the Majorana prob- 
lem. Depending on the parameters, the model can realize dif- 
ferent gapped and gapless phases. The latter generally have 
Fermi surfaces, and here we are focusing on such gapless 
phases and their qualitative properties. For all illustrations 
below, we use parameters {J x , J y , J z , J' x , J' y , J' z , h} 
< 1 .7. 1.4, 0.4, 1.0, 1.3, 0.2, 0.95} with x, y, z, x', y', and 
z' defined in Fig. Q] Gapless phases with Fermi surfaces ap- 
pear in wide parameter regimes, and we remark that there is 
no fine tuning of parameters to find such phases. The rea- 
son we choose to present the specific parameters is that in this 
case, the Fermi surfaces are sufficiently small, so when we an- 
alyze the singularities in the structure factors in Sec. IIIIBI it 
is easier to clearly see the locations of the singularities. 

For an illustration of how these two bands of usual com- 
plex fermion fields vary with momentum fe, we show them in 
Fig. |2(a)| along a cut with k y = — 37r/4. We label the bands 
from top to bottom as 1 to 2. We can see that only the band 2 
crosses the zero energy, which is true also when we scan the 



whole k y axis, and the populated Fermi pocket in the B.Z. is 
shown shaded in Fig. |2(b)[ ^i 

It is interesting to discuss qualitatively the thermodynamic 
properties in this phase. Because of the presence of the gap- 
less Fermi surface, such spin liquid is expected to show metal- 
like specific heat and spin susceptibility at low temperature, 
although the Wilson ratio is different from that of spin- 1/2 
fermions. 34 Furthermore, magnetic impurities coupled to this 
model would possibly show an unusual Kondo effect and 
Ruderman-Kittel-Kasuya-Yosida interaction^ 



III. CORRELATION FUNCTIONS 

We mainly focus on the spin correlations and spin-nematic 
correlations, and since there are four sites per unit cell, there 
are many correlation functions one can define. However, since 
all the spin correlations show similar behaviors among them- 
selves and so do spin-nematic correlations, we consider spe- 
cific examples defined as 



with 



^(r) = (S+(r,2)S-(0,2)), 
7- 2 (r) = <P+(r)p-(0)>, 



S +/ - = S X ± iS y = {<J X ± ia v )/2, 



(15) 
(16) 



(17) 



and magnon-pair creation operator 

P+(r) = 5+(r,2) 5+(r + x,4) 

= 5+(r c -|,2) 5+(r c + |,4), 



(18) 



associated with the bond (r, 2; r + x, 4). In the last line, 
r c = r + x/2 is the center of mass coordinate, and £ is the 
vector joining the two sites of the nematic operator, which 
is simply £ = x here. The magnon-pair operator, Pf, = 
S+S+ for a local pair of sites {j, k}, describes spin-nematic 
properties^— and can be connected to the usual traceless 
rank two quadrupolar tensor defined as 



1 



Sj ■ s k 



(19) 



through P+ = Q--Q. fc 

Furthermore, power-law correlations in real space corre- 
spond to singularities in momentum space, which we can 
study by considering the corresponding structure factors 

2? 1/2 (q) = ^J- 1/2 (r)e-^ r . (20) 



A. Long wavelength analysis 

We focus on the long-distance behavior and therefore retain 
only the contribution from band-2. The spin operator can be 
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(a)Bands of complex fermions, f^' y ' z (k), along a cut at k y = —3iz/4. 
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(b)Contour plot of band 2 of complex fermions, /if ,y,2 (k). 



FIG. 2. (a) Illustration of energy spectra of the two bands 
of the complex fermion fields along a cut with k y — 
— 37r/4. Here we take parameters {J x , J v , J z , J' x , J' y , J' z ,h} — 
{1.7, 1.4, 0.4, 1.0, 1.3, 0.2, 0.95}. Diagonalizing the Aij matrix in 
Eq. j 1 2b gives four bands, but only the two bands shown are needed 
for solving the Majorana problem (see text). Because only band 2 
crosses zero energy, we simply focus on it for long wavelength anal- 
ysis, (b) Contour plot of band 2, with the occupied Fermi pocket 
shaded. 



compactly written as 



S a (r,a)~ 

E E \N^{a)^y^(\,)p 2 ^)e-^-^ 

k.k'eB.Z. 8,-y ^ 



where Mkk' 



-W2,k(a)^2,k'(a)/iV MC and N\ 



kk' 



7 2,k 



{a)v 2X {a)/N uc = -N^. 



In order to determine long-distance behavior at separation 
r, we focus on patches near the Fermi surface of band 2 where 
the group velocity is parallel or antiparallel to the observation 
direction n = r/\r\, because at large separation |r| kp , 
the main contributions to the correlations come precisely from 
such patches. Specifically, we introduce Right(R) and Left(L) 
Fermi patch fields and the corresponding energies 



fp' ifl) (Sk) 



J2 K^FP 



■5k) 



(Sk) 



) FP\ 



P8h\ 



-6k± 



(21) 
(22) 



where the superscript (n) refers to the observation direction 
and P = R/L = +/— ; Vpp is the corresponding group ve- 
locity (parallel to h for the Right patch and anti-parallel for 
the Left patch); <tp = R/ l is the curvature of the Fermi surface 
at the Right/Left patch; Sku and Sk± are respectively compo- 
nents of 6k parallel and perpendicular to n. It is convenient 
to define fields in real space 



fp V) 



E 

<5fc£ Fermi Patch 



fp' {il \5k)e 



5k-- 



(23) 



which vary slowly on the scale of the lattice spacing [and 
from now on we will drop the superscript (n)]. In this long- 
wavelength analysis, the relevant terms in the spin operator 
are 

S a (r,a)~ 

~ E E {^PP'(«)^ 7 /p t ( r )/F'( r )e" i(kFP " kFP ' ) ' r + 

P.P' B.-y ^ 



+ 



P.P' B,"t 



(24) 



The above long-wavelength expression for the S a opera- 
tor implies that the corresponding correlation function defined 
in Eq. (fT~5T > contains contributions with q = 0, \tFR . — k F £, 
±2kp, and ±(kpij , + \s.fl)- More explicitly, for a patch spec- 
ified by ep(8k) in Eqs. (f2TT>-(l22T>. we can derive the Green's 
function for the continuum complex fermion fields as 

Using this and Eq. d24l . we can obtain the spin correlation 



•Fi(r) 



\Nrr\ 



\Nll\ 



(26) 
(27) 



2\N RL \ 2 sm[(k FR -'kFL) t] 

4 /2 4 /2 |r| 3 
\M RR \ 2 sm(2k FR -r) \M LL \ 2 sin(2k FL ■ r) 

£fl|r| 3 C^rf ( ' 

2|M flL | 2 cos[(k Ffl + k FL )-r] 

£r £l M 3 
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where we used N LR — —N* RL and M LR = M RL . 

For the long-wavelength description of the spin-nematic 
correlations, we can in principle plug the expression of spin 
operator, Eq. (l24l . into either Eq. £00 or Eq. ([19). We re- 
mark that even though the microscopic spin-nematic operators 
contain four local Majorana fermions expressed in general as 
c ? c k c ] c k> wnen calculating the correlation functions, there 
are cases when pairs of Majorana fermions Wick-contract lo- 
cally and produce a constant factor. Take Q xy as an example, 
Q 



( c j c k + c j c k) c j c k> an< ^ observe that the last two Ma- 
joranas c|c| can Wick-contract when calculating the correla- 
tion functions. For this reason, the spin-nematic correlations 
show the same dominant power-law behavior as spin correla- 
tions, and effectively we have fermion bilinear contributions 
to the spin-nematic. 

From now on, we focus on the dominant contributions to 
the spin-nematic correlations. The diagonal and off-diagonal 
elements of the quadrupolar tensor can be written concisely 
using center of mass and relative coordinates, where we define 
(j, k) = ({r, a}, {r' ; a'}) = ({r c - £/2, a}, {r c + £/2, a'}), 



s aa ~EE 

PP' P^a 



Af, a p,e l(kFP+k ^' )r <=/p/p' 



H.c.k (30) 



_1_ R°°' „-'»(l!FP-k F p/)-r c fP\ fP 
T-Dpp'6 Jp Jpi 

Qaf> ~ 111 { 4p'e ,(kFP+k ^' ),rc (/p/^+af>/3) 
PP' ^ - 



-B P a p,e- i(kFP " k ^' ) - rc (/p t /p' + a+*0) 



H.c. }.(31) 



aa 
PP' 



-iv 2 A a >2,P'(a')e-^r- k F P '>Z/y Nuc 
-w*.p(a)u 2 ,p,(aV (kFP+kFP ' K/ 7^« c ; the 



Here, A 

and B a p a p , 

slowly varying fermion fields are evaluated at r c . The above 
implies that the spin-nematic correlations contain dominant 
contributions at wavevectors q = 0, kpp — kp^, ±(kpp + 
kpi). Note that the contributions with q = ±2kp vanish by 
Fermi statistics. 

The long-wavelength expression for the dominant contribu- 
tions to the spin-nematic correlations, Eq. ([TBI , is 



^ 2 (r) 



(B 



RR 



Bt 



RR) 



\B RL + B 



2<£p|r c |3 

2 



(32) 



LP I 



■ sin[(kpp - kp L ) • r c ] (33) 



irl 



A 



LR. | 



cos[(kpp + kp L ) • r c ] , (34) 



where we abbreviate A 



a=2;a'=4 
PP' 



Appi and B P p 2 ,' a 4 



Bp pi. The above long- wavelength descriptions can be used 
to analyze the data obtained by exact numerical calculations. 
Here we also note that the model does not have time-reversal 
and inversion symmetries, so the location of the correspond- 
ing R-L patches which are parallel or antiparallel to the obser- 
vation direction can not be determined easily and need to be 
found numerically. 



Before leaving this subsection, we remark that we can sim- 
ilarly analyze local energy operators such as bond energy,^ 



J 3 J 3 - 



the long wavelength description con- 



tains terms fpfp, and Y!,p fp'fp' tnat ?ae spin-singlet 
variants of terms in Eq. (l30l . We can therefore see that the 
spin, spin-nematic, and local energy observables cover all 
fermionic bilinears. Below, we focus on the spin and spin- 
nematic operators. 



B. Exact numerical calculation 

We calculate the spin correlations, Eq. (Q3), and the spin- 
nematic correlations, Eq. ( fT6] l, for any real-space separations 
r and confirm that they have the same dominant power law 
envelope l/|r| 3 . For an illustration, we show the spin correla- 
tions and spin-nematic correlations for r along a specific di- 
rection, e.g. i-axis, calculated on a 300x300 lattice. In Fig. [3] 
the log-log plot of | J 7 ! (r) | and | Ti (r) | clearly shows the same 
1/ 1 r | 3 envelope. In addition, the irregular behavior of the data 
is due to oscillating components. The wavevectors of the real- 
space oscillations form some singular surfaces in the momen- 
tum space, which we analyze next. 

Focusing on the structure factors T>i/ 2 (q) defined in 
Eq. d20"l i, we calculate the spin correlation and spin-nematic 
correlation at each site within a 100x100 lattice and nu- 
merically take Fourier transform. Figure |4(a)| gives a three- 
dimensional (3D) view of the spin structure factor. We can 
clearly see cone-shaped singularity at q = 0, which is ex- 
pected from Eq. 



Z>i(q~0)~ |q| 



(35) 



A closer look at the spin structure factor also reveals sin- 
gular surfaces at kpp — kp^, ±2kp, and ±(kpp + kp^), 
as expected from Eqs. ([2Tb-(l29b. In order to see the loca- 
tions of the singular surfaces more clearly and compare with 
our long-wavelength analysis, we show top view of T>\ (q) in 
Fig. |4(b)| We numerically calculate the wavectors 2kpp and 
Q ± = kpp ± kpi for all observation directions (by first 
finding corresponding Right and Left Fermi points with anti- 
parallel group velocities) and superpose the traced lines on the 
figure. We see that the lines we get from the long-wavelength 
analysis match the singular features in the exact spin structure 
factor. Note that the singularities are expected to be one-sided, 

P 1 (Q_+«5q)~|(5g||| 3 / 2 e(-^||), (36) 
V^kpp +Sq) ~ |^||| 3 / 2 e(-5 g ||) , (37) 
2MQ + + <5q) ~ |<fy|| 3 / 2 e[-5g||Sign(£p -££.)]. (38) 

The first and second equations are singular from the inner side 
of the central "ring" and the closed rings sitting roughly on 
one diagonal of the B.Z. in Fig. |4(b)| and the last equation is 
singular from the inner side of the small "triangles". 

Similar analysis can be applied to the spin-nematic struc- 
ture factor except that there are no ±2kp singularities. The 
3D view of the spin-nematic structure factor is shown in 
Fig. |5(a)| and we can clearly see the q = singularity and 
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(a)Spin correlation, Eq. U51 . for r = xx 
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(b)Spin-nematic correlation, Eq. | |16K for r = xx 



FIG. 3. Figures (a) and (b) illustrate power-law behaviors of the 
spin correlations and spin-nematic correlations. We calculate J-i(r) 
and .7-2 (r) with r taken along the £-axis for a system containing 
300x300 unit cells. The log-log plots in (a) and (b) clearly show 
the same dominant l/x 3 envelope (straight line in the figures). Here, 
we show the absolute values of |J-"i(r)| and |J-2(r)|, and indicate 
the sign with open square boxes for negative correlations and filled 
circles for positive correlations. The irregular behaviors are due to 
oscillating parts. The reason that the panel (a) has no positive data 
is likely because the non-oscillating part is quantitatively stronger, 
although we can still see oscillations about the 1 /x 3 line. 



q = kpR — k.pL singular line (central ring). The q = 
±(kpfj + kpL) singular lines (small triangles) are quite weak 
but still visible, and their locations can be seen more clearly 
in the view from top shown in Fig. |5(b)| 



IV. MAJORANA SPIN LIQUID IN THE ZEEMAN FIELD 

In the presence of the Zeeman magnetic field, we need to 
consider the additional term in the Hamiltonian 

Hz = ~B z J2 a i = B z Y,^ x (r,a)c y {r,a), (39) 

i r,a 

where we used explicitly the rewriting of spin in terms of Ma- 
joranas, Eq. ([6j. We remark that the Zeeman magnetic field 




7T 

(a)Three-dimensional view of the spin structure factor 




(b)Top view of the spin structure factor 

FIG. 4. (a) 3D view of the spin structure factor, ©i(q), defined in 
Eq. i20\ . We can clearly see the singularity 2?i(q) ~ |q| at q = 
and we also see weak singular lines: one forming central closed 
ring; two closed rings sitting roughly on one diagonal of the B.Z.; 
and additional weak singular features near the centers of the latter 
rings, (b) These singular lines are brought out more clearly when the 
structure factor is viewed from top. We superposed the locations of 
the singularities calculated using the Fermi surface information: The 
inner blue ring specifies the line at ]<.fr ~\<.fl\ the red closed rings 
specify the lines at ±2kf; the small green triangles specify the lines 
at ±(k F R + k F i). 



only couples to the spin degrees of freedom (c a Majoranas) 
and not to the orbital degrees of freedom {d a Majoranas that 
produce the Z-2 gauge fields). Therefore, the model is exactly 
solvable even in the presence of the magnetic field. Through- 
out, we assume the K p terms, Eq. (Q]i, are large enough so that 
the ground state remains in the zero Z2 flux sector. 
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The Hamiltonian in the Zeeman magnetic field becomes 



(a)Three-dimensional view of the spin-nematic structure factor 



% 





(b)Top view of the spin-nematic structure factor 

FIG. 5. (a) 3D view of the spin-nematic structure factor, 2?2(q), 
denned in Eq. d20l l. Like the spin structure factor, there is a clear 
singularity X>2(q) ~ |q| at q = and there are weak singular lines, 
one forming central closed ring and the other forming small triangles 
sitting roughly on one diagonal of the B.Z.. (b) The view from the top 
shows more clearly the location of the singular lines. We superposed 
the locations of the singularities calculated using the Fermi surface 
information: The inner blue ring specifies the line at \lfr — ^fl, and 
the small green triangles specify the lines at ±(kp.R + \lfl)- Note 
that unlike the spin structure factor, there are no ±2ki? singularities. 



It is interesting to note that the Zeeman term only af- 
fects the c x and c y Majoranas while leaving the c 2 Majorana 
unalteredi^i We can diagonalize the Hamiltonian by start- 
ing with the zero-field solution, Eqs. (fl3l l- ([l4l) . We define two 
complex fermion fields, 



n = J2 E i^(k) + 2B z ) 

6=1 kGB.Z. 

2 

+ E E (2cb(k)-2B,) 



6=1 kGB.Z. 

E E 2 < b (k) 

6=1 kGB.Z. 



/t_(k)/ tl _(k) 
/J + (k)/ 6 ,+(k) 



/r(k)#(k)-- 



(41) 
(42) 
(43) 



Jb. 



b,+/ 



.(k)=[tft(k)±i#t(k)]/>/2. 



(40) 



This form implies that / T + carries S z quantum number +1 

and fl _ carries S z = — 1, while f z ^ carries S z = 0. 

An interesting property in this model is that the f z Fermi 
surface (associated with the c z Majorana) remains no matter 
how large the magnetic field is>2i therefore, there are always 
gapless excitations in this system. 

For an illustration of several different phases that can oc- 
cur under the magnetic field, we take the same parameters as 
in Fig. |2] and examine the effective Zeeman shifting of the 
band-1 and band-2 for each complex fermion species. Fig- 
ure [6] shows the magnetization as a function of magnetic field 
B z . There are quite rich features in this model. When we turn 
on the magnetic field, the /j, + bands move downwards while 
the _ bands move upwards. First, the field increases up to a 
threshold value, roughly B z = 0.024, where the Fermi surface 
of the fi- vanishes and there is a discontinuity in the slope of 
the magnetization curve shown in the inset in Fig. [6] The f% t + 
band is pushed down and the fy + Fermi sea keeps growing 
until it completely covers the B.Z. at B z ~ 0.5. For the field 
between 0.5 and 2.8, the f\ t+ band remains above the zero 
energy, and we have the half -polarized magnetization plateau 
phase. The band reaches zero energy at B z ~ 2.8, we 
leave the first plateau phase and the magnetization starts to in- 
crease. When the magnetic field is large enough to completely 
push the /i + band below zero, B z ~ 3.3, the /i + Fermi sur- 
face also vanishes and we enter the fully-polarized phase, the 
second plateau phase. We also see some weak features in the 
regimes of increasing magnetization that are due to the van 
Hove singularities when the energy passes the saddle points 
of bands 1 or 2 [cf. the contour plot of band 2 in Fig. |2(b)| , 
but these van Hove singularities are rather weak in 2D. 

We remark that in all regimes, the Fermi surface of the /| 
remains the same and gives gapless excitations. The fully po- 
larized phase is actually the original Kitaev-type model pro- 
posed by Baskaran et al. 41 This can be seen either directly by 
examining the physical Hamiltonian Eq. ([]]), or in the Majo- 
rana representation where <j z — 1 = —icjCj, so the constraint 
becomes D; = —ic^c^c^d^d^d 2 . = c z d x ,d v -d z = 1, and the 

J J J J J J J J J J J 

model in terms of the orbital degrees of freedom reduces to 
thatinRef. JH. 

It is interesting that in the half-polarized plateau phase, 
even though the spin excitations are gapped, the spin degrees 
of freedom are entangled in the ground state. The spinless 
gapless excitations can be in principle detected by measuring 
bond-energy correlations or by entanglement entropy calcu- 
lations. Finally, the regimes of increasing magnetization can 
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FIG. 6. The magnetization per site, M z = (^2 . o"J)/iV s it s, can 
be calculated as (pi,+ + p2,+ — pi,- — P2,-)/2, where pb,± is the 
population of the fb.± fermions per unit cell. The magnetization 
curve for our model in the Zeeman field shows rich features (see text) 
including half- and fully-polarized plateaus. The inset shows a blow- 
up of the small B z region where the Fermi surface /2,_ disappears 
and there is a slope discontinuity. We note that the magnetization 
inside the unit cell is not completely uniform (because of the reduced 
lattice symmetries) but is quantitatively similar at each site. 

be viewed as generic compressible Bose-metals 8 in the model 
with global U(l) symmetry. 



venient tool to reduce the problem to more familiar calcula- 
tions. However, if the symmetries of the model are sufficiently 
low, there can be an "emergent" conservation of f-s. For ex- 
ample, if in our treatment we have a very small Fermi pocket 
of f-s centered at some momentum K = (K x ,K y ), then the 
four-fermion terms /'/ / / carry approximate momentum 
4K while / V W carry « 2K, so these terms are not allowed 
by momentum conservation if K is some generic non-special 
wavevector. In this case, only the familiar terms PPff con- 
serving the fermion number are left and we have the emergent 
fermion number conservation law in the long-wavelength the- 
ory. The above is also true in our paper in which the non- 
symmetric small Fermi pockets are realized. Furthermore, 
only benign forward-scattering four-fermion terms survive in 
our model since the Cooper-pair interactions carry non-zero 
momentum. We do not need to consider six- and eight- 
fermion terms as they are irrelevant in the Renormalization 
Group(RG) sense. 

We have chosen a model that lacks time-reversal and lat- 
tice point group symmetries so as not to worry about possible 
residual pairing instabilities away from the exactly solvable 
limit. Such instabilities can be relatively weak also in more 
symmetric models and the discussed phenomenology can ap- 
ply in these cases as well. It should be possible to explore the 
stability and the nearby phases by studying such models on 
ladders^ using weak-coupling RG technique and Bosoniza- 
tion analysis^— 



V. DISCUSSION 

We proposed SU(2)-invariant Kitaev-type model on the 
decorated square lattice that realizes QSL with parton Fermi 
surfaces. Having the benefit of the exact solutions, we can 
make robust general observations about such SU(2)-invariant 
Majorana QSL. One of the distinguishing characteristics of 
this state is that it has strong spin and spin-nematic fluctua- 
tions as manifested by the same power-law behavior in the 
correlations considered in our work. Because of the finite 
density of states at the Fermi surface, properties such as the 
specific heat, spin susceptibility, and NMR relaxation rates 
are essentially similar to a metal. Note that the volume/shape 
of the Fermi sea can be arbitrary, and we can easily tune the 
model to have larger or smaller Fermi pockets. In this way, the 
phase is quite distinct from the conventional spinon Fermi sea 
QSL^i where half of the Brilloin zone is populated by spin- 1/2 
spinons before the Gutzwiller projection. Given the variabil- 
ity of the Fermi surface, the properties of the Majorana QSL 
can be very sensitive to parameters, which we can also easily 
tune to produce gapped phases. 

Let us briefly discuss stability of the exactly solvable model 
to general perturbations. First, we note that our complex 
fermions f a are not conserved microscopically. In principle, 
allowed four-fermion interactions would contain terms such 
as fafpfifl and flfpfjfS' in addition to the more familiar 
terms f^f^fjfs- The origin of the non-conservation of the 
complex fermions is because the microscopic fields are Majo- 
rana fermions, while we used the complex fermions as a con- 



In the presence of the Zeeman magnetic field, there are 
more interesting phases with distinct stable Fermi pockets of 
fb,+t fb,-> and f§ ■ The Zeeman field breaks the global SU(2) 
down to U(l), and the compressible phases in the field are 
Bose-Metal-Mke phases&S 

We also found an interesting plateau phase at half- 
magnetization. Due to the gap for spin excitations, it is a spin 
insulator, but since the f z remains gapless, we expect to still 
have metal-like specific heat and thermal conductivity. 

This behavior in the plateau phase arises because some of 
the parton constitutents of the spin operator acquire a gap 
(band gap in the present case). Some such physics per- 
haps can be relevant for the explanation of very recent NMR 
experiments 24 in EtMe3Sb[Pd(dmit)2]2 showing a drastic re- 
duction in the spin relaxation below temperature of the order 
IK as if a spin gap opens up, while the thermal conductiv- 
ity measurements and thermodynamic measurement s 25,27 are 
consistent with the presence of a Fermi surface of fermionic 
excitations down to the lowest temperatures. This phe- 
nomenology is also qualitatively similar to our recent paper. 7 -^ 
working in a setting closer to the EtMe3Sb[Pd(dmit)2]2 exper- 
iments. We considered a scenario in which, upon writing the 
spin operator as S + = f^~f±, there could be a phase in which 
one spinon species becomes gapped due to pairing, while the 
other species retains the Fermi surface. It is fascinating to 
further explore such idea where some partons are gapped and 
some are gapless in more realistic settings. 
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